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I. Abstract 

Wigner functions play a central role in the phase space formulation of quan- 
tum mechanics. Although closely related to classical Liouville densities, 
Wigner functions are not positive definite and may take negative values on 
subregions of phase space. We investigate the accumulation of these negative 
values by studying bounds on the integral of an arbitrary Wigner function 
over noncompact subregions of the phase plane with hyperbolic boundaries. 
We show using symmetry techniques that this problem reduces to computing 
the bounds on the spectrum associated with an exactly-solvable eigenvalue 
problem and that the bounds differ from those on classical Liouville distri- 
butions. In particular, we show that the total "quasiprobability" on such a 
region can be greater than 1 or less than zero. 



II. INTRODUCTION 

Since its introduction—, the Wigner function has been the subject of exten- 
sive study in the fields of quantum physics, quantum chemistry and signal 
analysis (se o 2 i 3 ' 4 i 5 ' 6 i 7 i 8 ' 9 i 10 and references therein). Since Wigner functions 
represent quantum states on phase space, they play a key role in the phase 
space formulation of quantum mechanics. They are also designed to closely 
resemble the joint densities of position and momentum, known as Liouville 
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densities, that are used in classical mechanics. In quantum physics, such 
studies have been stimulated in recent times by the development of quantum 
tomography, which has enabled the reconstruction of Wigner functions cor- 
responding to states of a variety of quantum systems^. Such experimental 
observations have confirmed that Wigner functions can be negative on sub- 
regions of phase space. This is one of several properties that can be used to 
distinguish Wigner functions from classical Liouville densities. 

The study of these "quantum properties" has been approached in a num- 
ber of ways including calculations of pointwise bounds on Wigner functions 
and bounds on various moment s 12 i 13 i 14 i 15 i 16 . A more recent development has 
been the study of bounds on integrals of the Wigner function over subregions 
of the phase spac o 17 i 18 i 19 , which we denote by T . We call such integrals 
quasiprobability integrals (qpis). For a given subregion S of T, the problem 
of determining best possible upper and lower bounds on all possible qpi's 
over 5* has been shown to be equivalent to the problem of determining the 
supremum and infimum of the spectrum of the region operator associated 
with S. This operator is just the image under Weyl's quantization map^i 
of the characteristic function of S, namely the function that equals 1 on S 
and elsewhere on T. In the special case of a quantum system with one 
linear degree of freedom, it has been shown that for any subregion of the 
phase plane enclosed by an ellipse, the eigenvalue problem is exactly solvable 
and the bounds on qpi's can be obtained analytically for ellipses of arbitrary 

17 

size*! 

The determination of bounds on qpis is important not only because it 
provides information about the structure of theoretically possible Wigner 
functions, which is a question of mathematical interest, but also because 
an understanding of that structure provides checks on experimentally deter- 
mined Wigner functions. It is therefore of interest to know if there are other 
subregions of the phase plane, and more generally of phase space, for which 
the spectrum of the associated region operators, and hence the best possible 
upper and lower bounds on all possible associated qpis, can be determined 
exactly. In this paper, we show that an exact formula for the spectrum of 
the region operator, from which the bounds are easily obtained numerically, 
can be derived for subregions of the phase plane with hyperbolic symmetry. 
The solvability of the eigenvalue problems for the corresponding region op- 
erators, as in the case of elliptical subregions discussed earlier, relies on the 
invariance of these regions under one-parameter subgroups of the metaplectic 
group Mp(2,M) of transformations of the phase plane. This group consists 
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of all real transformations of the form 

T : (q,p) -> (q',p) = (aq + /3p + g ,7<? + $P + Po) , (!) 

where a5 — /3'j = 1. In this paper, the subgroup of Mp(2,M) formed by the 
transformations T a : (q,p) —> (aq,p/a) , cr > is of particular importance. 

Several illustrative examples of eigenvalue problems for hyperbolic regions 
are considered in what follows, including the interesting limiting case of an 
infinite wedge. We shall be concerned with quantum systems with one linear 
degree of freedom, described in terms of a Hilbert space of states 7i, and 
with the properties of Wigner functions on the associated (q, p) phase plane 
T. We are not concerned with dynamics, and consider each Wigner function 
at a fixed time. Dimensionless phase plane coordinates (q,p) are used, and 
in effect we set h — 1. Finally, we note that in the absence of limits of 
integration, integrals are assumed to run from — oo to oo. 



III. BOUNDS ON QUASIPROB ABILITY IN- 
TEGRALS 

The Wigner function corresponding to a pure state ip E7i has the definition 

W^q,p) = ±-J^(q + r)i>(q - r)e 2 ^dr. (2) 

For a mixed state, the Wigner function is a convex linear combination of 
such integrals. It is known that Wigner functions are bounded at every 
point (q,p) E T such that —1/tx < W(q,p) < 1/n and that they satisfy the 
normalization conditions 



Wdqdp = 1 , < / W 2 dqdp < — 
Jr 



2tt 



where the value 1/2% is attained if and only if W corresponds to a pure state. 

More generally, an operator A is unitarily related to a phase space func- 
tion A(q,p) by the Weyl-Wigner transform— and its inverse, 

A = W(A), A = W-\A). (3) 

Here W -1 is Weyl's quantization map and W is such that the Wigner function 
corresponding to a quantum density operator p is given by W p = W(p)/(2%). 
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In this paper we make extensive use of the configuration realization, in which 
A can be expressed as an integral operator 

(Aip)(x) = J A K {x,y)^{y)dy . (4) 

We refer to the function Ax{x,y) as the configuration kernel of A. It is 
related to the phase space function A(q,p) by the formula s 22 ! 23 

A(q,p) = J A K (q-y/2,q + y/2)e ipy dy, (5) 

A K (x,y) = ^jA((x + y)/2,p)e i ^dp, (6) 

which provide an explicit realization of the transformations (J3J). 

An important property of Wigner functions is that quantum averages on 
phase space take the same form as classical averages: if A(q,p) is the phase 
space representation of a quantum observable A, then its quantum average 
in the state with density operator p and corresponding Wigner function W p 
is given by 

(A) = Tr(Ap) = J^W p (q J p)A(q J p)dqdp. (7) 

The qpi of a Wigner function W over a subregion S of T may be written 
as the functional 

Qs[W] = I W(q,p)dqdp. (8) 
Js 

Note that the integral on the RHS can be rewritten in terms of the charac- 
teristic function Xs{l,p) that equals 1 on 5 1 and on its complement: 



Qs[W] = J W(q,p)xs(q,p)dqdp, 

and, by comparing with (JZj), we can write 

Qs[W}=(xs), (9) 

where we have introduced the region operator xs — W^ 1 (xs)^S with config- 
uration kernel (as given by (JHJ)) 

XsA^v) = ^ / Xs((x + y)/2,p)e l ^-^dp. (10) 



Since the expectation value of a quantum operator always lies between 
the extremal values of its spectrum, we deduce from Q that Q,s[W] must lie 
between the infimum and the supremum of the spectrum of xs- Moreover, 
as the spectral bounds on the expectation value of an operator can be ap- 
proached arbitrarily closely with normalized states in 7i, these bounds are 
best-possible. Hence the best-possible bounds on the qpi functional Q$ are 
provided by the extremal solutions to the integral equation 

{Xsip){x) = J Xs,K(x,y)4>(y)dy = \ij){x) (11) 

that defines the eigenvalue problem for xs- 

For a general region S, the integral equation (jllj) is not exactly solvable 
and the bounds on its spectrum must be obtained by using computational 
methods. However, there is a subclass of regions for which the (generalized) 
eigenvalues and eigenfunctions can be determined exactly. This subclass is 
the set of regions that are each invariant under a one-parameter subgroup 
of the metaplectic (or linear canonical) group of transformations of the 
phase plane. Any such transformation U has the special property that its in- 
verse Weyl-Wigner transform U = W _1 (£7) is a unitary (and thus spectrum 
preserving) operator acting on 7i. If a subregion S of V is invariant under 
a metaplectic transformation, then it follows that the associated region op- 
erator xs is invariant under the corresponding unitary operation U, that is 
generated by an operator f of no greater than the second degree in q and p. 
It follows that [xs, r] = and hence that the eigenfunctions of xs ma y be 
chosen so that they are also eigenfunctions of f. These are readily obtained 
by solving the eigenvalue problem for f. 

This approach can be applied to regions that are bounded by ellipses, hy- 
perbolas, parabolas and straight lines. (If the boundary is composed of sev- 
eral curves, then each curve must be invariant under the same one-parameter 
subgroup of Mp(2,R).) In the case of elliptical regions, the best possible 
bounds have already been described—, while the fact that the marginal dis- 
tributions of the Wigner function are true probability density functions^ 
implies that integrals over regions bounded by parallel straight lines must lie 
in the interval [0, 1]. In this paper, we consider the problem of determining 
the best possible bounds on qpis over regions with hyperbolic boundaries. 
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(a) (b) 




Figure 1: Graphs of hyperbolic regions and their boundaries: in (a), the 
hyperbolic region Sk is shown as are the boundary curves and , while 
in (b), the infinite wedge So is depicted. 

IV. BEST POSSIBLE BOUNDS ON QPIS 
FOR HYPERBOLIC REGIONS 

In order to demonstrate our technique for constructing the bounds on qpis, 
we begin with a simple example. Let Ck be the hyperbolic curve consisting 
of all points that satisfy 

qp = k, k > 0, (12) 

as depicted on in part (a) of Figure [T] Note that Ck is itself composed of 
two curves, namely C£, which lies in the positive (q,p) quadrant of T and 
CjjT, which lies in the negative (q,p) quadrant of T. It is clear that the 
curves are separately invariant under the action of the transformation 
T a : (q,p) — »■ (aq,p/a) for all a > 0. 

Now consider the subregion Sk that contains all points in T such that 
qp > k, q > 0, which is indicated by the shaded region in part (a) of Figure 
[0 Since this region can be viewed as the union of all hyperbolic curves 
with I > k, it is itself invariant under the action of T a . In order to apply 
this symmetry to the problem of determining the bounds on qpis over Sk, we 
must first construct the corresponding region operator, which we denote by 
Xk- Note that the characteristic function on Sk, may be written as 

** = {l W oL™ S e° -*/*), (13) 
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where H is the Heaviside function. Using JBJ), the configuration kernel for 
the region operator can be determined (see the Appendix for details): 



' x + y 



\5{x-y) 



2ni(x — y) 



(14) 



and hence the bounds on qpis are given by the spectral bounds associated 
with the integral equation 



g x+y 



\8{x-y) 



1 



2iri(x — y) 



^{y)dy = fiijj(x) 



(15) 



We know, however, that the region operator \k is invariant under the 
set of operator transformations that correspond to the subgroup of Mp(2, M.) 
formed by T a , a > 0. Since the effect of T a is to squeeze position and stretch 
momentum (or vice versa) while preserving the canonical commutation re- 
lations, the corresponding operator transformation, up to an unimportant 
phase, is given by the squeezing operator U a = exp {i log a (qp + pq)/ '2). This 
implies that Xk commutes with U a for all a > 0, and hence 

[Xfc,w] =0, u = (qp + pq)/2. 

It follows that the eigenfunctions of \k can be chosen such that they are 
also eigenfunctions of Cj. We can then obtain a partial solution to the inte- 
gral equation (|15J) by solving the equation uip = uoip. A number of results 
connected with this problem can be found in a paper of Chruscinski— . On 
configuration space, this equation appears as the first order differential equa- 
tion 

(iw-ity. (16) 



x 



dx 



The solutions of this equation are complex- valued linear combinations of the 
functions 



\ e iulog|i| 

73^ M 1/2 




x > 
x < 



X > 0, 

iui log | x (17 ) 

x < . v ' 



V2i |x|V2 



Here u can take any real value. These solutions are generalized functions 
and are elements of the space of tempered distributions Q — , of which Ti is 
a proper subspace. The factor l/y/2n is inserted to ensure that (ipt^w 1 ) = 
5{uj — uj'). Since they have disjoint support, and ip~, are orthogonal for all 
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u),u' G R. Note that, since log \x\ — > — oo as |x| — > 0, the eigenfunctions be- 
come highly oscillatory in the neighbourhood of the origin and are undefined 
at |x| = 0, due to the \x\ 1 / 2 term in the denominator. 

Since the ip^ form two independent families of solutions to (|16|L the 
eigenfunctions of Xk are not yet fully determined. In order to construct these 
solutions, we must solve the reduced eigenvalue problem 

XfcV'a, = fi(w,k)l() u , Vw = Ow^ + A;^, (18) 

where a w ,P u G C. In order to solve (fTHj) . we must first determine the action 
of Xfc on the two-dimensional subspace Q' u of £?' spanned by ip^ and ^J: 
Xki^t ^Z) T = A(u, ^){'4 > tji' l l ) Z) T i where A(ou, k) is given by the matrix 

V. A 2 i(^,fc) A 22 (w, fc) / 

The matrix elements of A(u, k) can be computed by using the configura- 
tion realization of Xk, details of which are presented in the Appendix. It so 
happens that the matrix elements of A(u, k) depend on the functions d(u, k) 
and a(u>, k), that are given by 



d(u,k) = - 



f r giioz— 2ik coth(«/2) 

tanh(7ra;) H 3" < (b — — — — dz 



An i J c cosh(2;/2) 



(20) 



where Co is any closed path in the complex plane that contains only the pole 
at z = 0, and 

e~™ / r / e ^-2fetan(t/2) cos(wt - 2A; tanh(t/2))\ 
~2T~ v7o V sin(t/2) sinh(t/2) / * 

°°cos(^-2fctanh(t/2)) 



sinh(t/2) 

The above formula is written in this way, because the individual terms in the 
first integral are singular at t — 0, whereas their difference is well-defined. In 
terms of these functions, we may expand A(lj, k) as 



A(w,k) 



\ + d(cu, k) \ [a(u, k) + ie~™{\ + d{u, k))] 

\[a{u,k)-ie-* w (l + d{u;,k))] 

(22) 
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Figure 2: Graph of the best-possible bounds on qpis over Sk for k in the 
range to 5. The inset graph, with k in the range to 0.2, shows that the 
upper bound lies above 1, but converges rapidly to 1 as k increases. 



Hence the spectrum of the region operator Xk splits into positive and negative 
parts, which we label by /i+(uj, k) and k) respectively: 



Of particular interest are the functions = ini ue nL[/,-(u,k) and Uk = 
sup w6R /i+(co>, k), since they provide the best-possible bounds on qpis over 
the hyperbolic regions Sk- Although it does not seem possible to obtain ex- 
act expressions for these functions, it is not difficult to compute the bounds 
after first evaluating a and d numerically. 

These bounds are graphed in Figure 121 for k in the range [0, 5] from which 
we conclude that the upper bound on qpis over Sk remains close to but 
greater to 1 for all k and that this difference is greatest when k = (see 
inset). The lower bound displays a more marked difference from the classical 
bound of 0, reaching a minimum value of —0.3089 at approx. k = 1.9, before 
rising again. A surprising result is that the lower bound does not approach 
for large values of k. Nonetheless, this appears to be a characteristic feature 
of bounds on qpis for many classes of regions^. 



1 1 



+ d(u, k) ± + d{u, /c)) 2 (l + e- 2 ™) + a{u, k) 2 



2 2 



(23) 



9 



IV.. 1 THE INFINITE WEDGE 



An interesting subclass of hyperbolic regions is provided by taking the limit 
as k — > 0. The region So so obtained is precisely the positive (q,p) quadrant 
of r, as depicted in part (b) of Figure^ Note that when k = 0, the functions 
d and a take a simplified form: 

d(u,0)) = itanh(7ru;), a(u,0) = — \ tanh(7Ro>) + u{uj) , (24) 

where u(u) may be expressed as an infinite sum^: 



u(u) = — > 



7T 



a -(uj 2 + (An + 1) 

n.=0 v ' 



If we now apply these simplifications to the spectral formula (I23j) . we obtain 
the spectrum for x : 

H±{u, 0) = | ^1 + tanh(vrcj) ± y/(2u(u) - tanh(vrcj)) 2 + (1 + tanh(vrc«;)) 2 ) , 

(25) 

which is graphed in Figure El The infimum and supremum can then be 
determined numerically, and to an accuracy of ±5 x 1CT 10 , we have that 

-0.155939843 < Q [W] < 1.007679970. (26) 

An interesting point is that these bounds are also best-possible when 
applied to regions defined by infinite wedges. This equivalence is due to two 
factors: firstly, by an appropriate metaplectic transformation T, the region 
So can be transformed into any infinite wedge with half- angle a < ir/2. 
Secondly, the operator transformation Ut that corresponds to T is unitary, 
and thus the spectrum of Xo is preserved under its action. This implies that 
the spectrum of any region operator corresponding to an infinite wedge is 
given by (|25|) . As a consequence, the integral of a Wigner function over any 
infinite wedge must lie between the bounds given in (|26|). Bounds on the 
spectrum of similar operators have been considered before, in the context of 
the quantum phase operator— and in connection with studies of probability 
backflow— but not to the same level of precision. 
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Figure 3: Graph of the spectrum of Xo which, as labeled, splits into the 
curves /i + (cj,0) and /i_(u;,0) 



V. EXAMPLES INVOLVING TWO BOUND- 
ARY CURVES 



As a second example, consider the slightly more complicated case of a region 
with a boundary composed of two curves with T a symmetry. There are several 
possible forms that such a region can take^, however, we will concentrate on 
just the subcase for which the boundary curves lie in positive and negative 
(q,p) quadrants, as shown in part (a) of Figure EJ 

In order to further simplify matters, we assume that both curves are 
labeled by the variable k. We label the class of regions that remain by Sk,2 
and note that this region may be written in terms of the region Sk as 



where R w denotes a rotation through an angle ir. Note that R n : (q,p) — > 
(—q, —p) and that the operator that corresponds to this transformation under 
the Weyl-Wigner transform is just the parity operator P, which acts on the 
canonical coordinate and momentum operators according to 



Sk,2 — Sk + R 7r (Sk) , 



(27) 



PqP 



q , PpP = —p . 
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Figure 4: Regions bounded by two hyperbolic curves: in (a), the hyperbolic 
region Sk,2 is shown while in (b) the double wedge 5*0,2 is depicted. 



Due to the linearity of the Weyl quantization map, this implies that the 
region operator that corresponds to 5^2 may be expressed as 

Xk,2 = Xk + PXkP ■ (28) 

This operator also commutes with uj, and hence its eigenstates can be chosen 
such that they are some linear combination of ip^. In order to find the correct 
linear combination, we must first determine the matrix representation of Xk,2 
on the subspace spanned by ip^. This turns out to quite simple, since the 
action of P on this subspace is given by Pip^ = VJ- Thus the matrix 
representation of Xk,2 is given by 

(29) 

The simple form of this matrix representation leads to the following ex- 
pression for the spectrum of Xk,2- 

A*2,±(w, k) = - + d(u, k) ± \a(u, k)\ , (30) 

In this case, the eigenf unctions are odd and even combinations of and 
are independent of k: 

^(x; k) = iifjx) = ±= (^+(x) ± ^p-(x)) , (31) 

which indicates that the operators Xk,2 commute for all k > 0. 
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Figure 5: Graph of the best-possible bounds on qpis over S^ 2 for k in the 
range to 5. 

The properties of the spectrum in this case vary somewhat from the 
preceding example. In particular, fi 2 _(ou,k) is not restricted to negative 
values and, similarly, /x 2 ,+ (u;, k) is not strictly positive, although clearly the 
inequality fj, 2 ,+(u, k) > fi 2 ~(uj,k) holds for all u G M, k > 0. Since the 
bounds on qpis over Sk,2 are given by the infimum Lk,2 and supremum U^ 2 
of the spectrum of Xk,2, it is these functions that are of primary importance 
in the context of this paper. Again, closed-form expressions do not appear 
to exist, so we must resort to computational techniques in evaluating these 
functions, the results of which are graphed in FigureEl In this case, the upper 
bound is well in excess of 1 for small values of k, but rapidly approaches 1 
as k increases. The lower bound, on the other hand, dips initially, reaching a 
minimum of —0.4014 at k = 0.4 before rising again, and appears to approach 
a finite negative value near to —0.3 as k — > oo. 

V..1 DOUBLE WEDGES 

It is again of interest to consider in more detail the limit as k — > of the 
region Sk t2 . The resulting region Sq >2 is the union of the positive and negative 
(q, p) quadrants (as shown in part (b) of Figure El and one might guess that 
qpis over such a region should be positive^, since it appears to composed 
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from the union of a set of infinite straight lines, over which the integral of the 
Wigner function is known to be positive. However, since these lines cross at 
the origin one cannot immediately apply this result and it will be shown that 
the true bounds on qpis lie significantly outside the [0, 1] interval to which 
classical probabilities are restricted. 

The region operator that corresponds to x.0,2 m ay be expressed in terms 
of x as Xo,2 = Xo + PXqP- The spectrum for this operator can be derived 
from (JBTij) upon substitution of (|24j) . from which we obtain 

tanh(7ru;) 



1 + tanh(7rc<j) 

02,±(",O) = J^^ ± 



u{uS) 



(32) 



This spectrum is graphed in FigurelHl and from this one sees that the function 
/i 2 = \ + u(oj) passes through both the infimum and the supremum of the 
spectrum of Xo,2- Since u(u) is an odd function, we need only calculate its 
global maximum in order to determine the bounds on qpis over 5*0,2 ■ Using 
computational techniques, this value can be obtained to great accuracy, and 
we find that the best-possible bounds (accurate to ±5 x 1CT 10 ) on qpis over 
5*0,2 are 

-0.236823652 < Q a , 2 [W] < 1.236823652. (33) 

Note that the upper and lower bounds sum to 1 since they are symmetric 
about 1/2. This symmetry can be explained by noting that if one rotates 
the region 5o,2 through an angle tt/2, then one obtains its complement: i.e. 
-^71-/2(50,2) = T/5o,2 = 5q 2 . Note that the integral of a Wigner function 
over r = S U 5 C is normalized to 1. Now, since the operator equivalent of a 
rotation is a unitary transformation, the region operator that corresponds to 
the complement of 5o,2 has precisely the spectrum given in (|32|). Accordingly, 
the spectrum of Xo,2 must consist of pairs that sum to 1 and, in particular, 
the upper and lower bounds on this spectrum must be symmetric about 1/2. 

As in the case of the region S , the bounds on qpis over 5 0> 2 can be ap- 
plied to a wider class of regions. We shall refer to elements of this wider 
class as double wedges, since they are formed by taking the union of an 
infinite wedge with its rotation through an angle it. By applying the ap- 
propriate metaplectic transformation, we can transform Sq^ into any double 
wedge. The corresponding operator transformation is unitary and preserves 
the spectrum of Xo,2, so that the spectrum of any region operator corre- 
sponding to a double wedge is given by (|32[). Accordingly, the integral of any 
Wigner function over an arbitrary double wedge must satisfy the inequality 
given in (|33|) . 
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Figure 6: Spectrum of the double wedge operator xo,2- 

VI. CONCLUSION 

The problem of constructing best possible bounds on integrals of the 
Wigner function is not only of mathematical interest, but should be of prac- 
tical significance in providing checks on experimentally reconstructed quan- 
tum states. Since our approach to the problem relies on specifying the region 
to be integrated over, it is important to identify the types of region for which 
the bounds can be easily computed. In this paper, we have considered several 
examples of regions with a hyperbolic symmetry for which the bounds can 
be computed numerically from the spectrum of an exactly solvable integral 
equation. We have demonstrated that the bounds on integrals of the Wigner 
function for these regions are not equivalent to those on integrals of true 
probability density functions. In particular, the lower bound is significantly 
below zero in all cases, although it lacks the scalloped effect arising from 
eigenvalue crossings as seen in the bounds for elliptical discs^. The upper 
bound also rises above 1 although for the most part the difference between 
its value and the classical bound is very small. This contrasts with the case 
of the disc, for which the upper bound always remains below 1. 

The results herein can also be extended to more complicated regions 
with boundaries given by an arbitrary number of hyperbolic curves sharing 
the same symmetry^, for example, the regions shown in Figure [H The 
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Figure 7: Examples of generalized hyperbolic regions are shown in (a)-(c). 

problem of determining the spectrum is essential the same but the matrix 
representations for operators corresponding to regions with many boundaries 
are functions of many variables and hence the behaviour of the bounds is 
much more difficult to characterize. 
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APPENDIX 

The configuration kernels that correspond to (|13|) take the form 

T7( E±y \ roc 

XkA^y) = A^ 1 J H(p-^)e^-^dp. (A-l) 

This integral can be computed in a generalized sensed, and we find that 

1 



Xk 



>K (x,y) = H(^)e 2ik ^ 



\5{x - y) 



2iri(x — y) 



(A-2) 



This expression for the configuration kernel of Xk enables us to determine 
the action of Xk on the space Q' u (recall that this is given by the matrix 
A(u, k) defined in (JUJ)) . In this representation Xk acts on as 



(xM(x) 



2ik^ 



\5{x-y) 



2ni(x — y) 



Mv)dy- (A-3) 



If we substitute = a^ip^ + (3^^ , then we find that the action of Xk when 
x < differs from that when x > 0. Thus, for x > 0, 



'o 

r-o 

+ Pu, e *+y 



and for x < 0, 



\5{x -y)- 
\5(x-y) - 



2m(x — y) 
1 

2ni(x — y) 



Juj log y 



e iw \og\y\ 

V^\y\ 



dy 
dy 



2ik^- 



\8{x - y) 



2iri(x — y) 



D iu) log y 



V^y 



dy. 



It is immediately clear that y^f^, k) = 0, since the x > case involves only 



The other integrals can be simplified and this process leads to the follow- 
ing expression for the matrix elements of A: 



A(u, k) 



~ + d(cu,k) \[a{uj,k) + ib{uj,k)\ 
| [a (a;, k) — ib(uj, k)} 



(A-4) 
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where the functions d, a and b are given by 

„ IN 1 f°° sin(u;t-2fctanh(t/2)) , 

d(u,k) = — / — t-, i \ JJ dt, A-5 

v ' ; 2tt J sinh(t/2) ' v ; 

, IN 1 f°° sin(cut-2fccoth(t/2)) , . . , 

2n J q cosh(r/2) 

,/ ,n 1 /"°° cos(c;t-2fccoth(t/2)) , 

Mw,fc = — / — i —. j- v 1 JJ dt. A-7 

v ' ' 2vr 7 cosh(i/2) v ; 

Note, however, that although the integrands of a and b are bounded for all 
t, they become highly oscillatory in the neighbourhood of the origin, which 
poses difficulties for numerical schemes. These problems can be alleviated 
by using the technique of contour integration. 

In the case of b(u,k), we consider the following contour integral in the 
complex plane 

i r p iuiz—2ikcoth(z/2) 

I c = 7T<P U /o\ dz . ( A - 8 ) 



2ir J c cosh(2;/2) 

where C is the contour shown in part (a) of Figure |H1 Although the contour is 
divided into four parts, only the integrals along the real axis contribute, since 
the contributions from the semi-circular segments vanish in the respective 
limits as e — > and R — > oo. Thus one has that 

poo giajt— 2ik coth(t/2) 

lim I c = 2— / — — - — dt, (A-9) 

a-oo ie -»o 2tt 7 cosh(t/2) ' v ; 

and result, 6(w, fc) = 3?{/c}/2. 

We can make use of the residue theorem in evaluating Ic'- 

I c = 2m Res/(«) , (A-10) 

where f(z) is the integrand in (|A-8|) . Note that f(z) has two distinct classes 
of residues: simple poles at z — (2n + l)ni and essential singularities at 
z = 2miri, with m G Z. The contour C encloses only the simple poles with 
n > and the essential singularities with m > 1, and hence the sum in 
flA-lOj) is over the residues at these points. 

It is easy to evaluate the residues at the simple poles and we find that 
the total contribution from the simple poles inside C is given by 

ReSs = J_^(_i)» e -aw-i)«w = sech ^) , (A-ii) 

n=0 
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(a) (b) (c) 




Figure 8: The contours used in evaluating the functions b(ui, k), a(cu, k) and 
d(uj,k): in (a), the semi-circular contour C, in (b), the contour C used to 
relate b(u, k) to d(u, k) and in (c), an example of a contour Co for evaluating 
R(u, k). 



The sum of the residues associated with the essential singularities can also 
be simplified: 

Res e = ±- V(-l)"e- 2 — k)== ^p^ R{u, k) , (A-12) 

n=l 

where R(ou, k) is the residue of f(z) associated with the essential singularity 
at the origin. 

Ordinarily, one might try to evaluate this residue by constructing the 
Laurent series for f(z), however in this case the coth term in the exponential 
makes this extremely difficult. However, one may use the residue theorem in 
reverse and evaluate R(u, k) by considering the integral of f(z) over a closed 
contour enclosing the origin (and no other poles): 

i r p iu}z—2ikcoth(z/2) 

R(u,k) = — <b , ; ■ dz. (A-13) 

v ' ; 2m J Cl cosh(^/2) v ; 

The rapid oscillations due to the coth term do not appear in this calculation, 
and due to its finite range this integral can be rapidly evaluated to a high 
degree of accuracy using numerical techniques. 

If we now collect the results for the residues together, we discover that 

sech(7ra;) 



b(u>, k) 



1 + — -%{R(u,k)} 



(A-14) 



Note that in deriving these results it has been assumed that uj > 0. A 
similar procedure (with the semi-circular contour defined in the lower half 
plane) enables us to extend the validity of (jA-14j) to all uj g R. 

1 
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In order to obtain superior expressions for the functions a(u, k) and 
d(u,k), we choose another contour (see part (b) of Figure [8J), this time 
involving five curves. However, we know from the above calculation that 
the integral over eo vanishes, which leaves four curves to consider. Of these, 
the integral over the imaginary axis from to n results in a pure imaginary 
contribution I n , while the integral over t\ contributes —1/2 in the limit as 
€\ — > 0. The integral over the positive real axis is equal to b(u, k) + a(u, k)i 
in the limit as eo — > 0, while the contribution from the line (ni, iri + oo) can 
be expressed as — exp(— ircu)(d(uj, k) + loot). 

By equating the real parts, we find that 

b(u, k) = e~™(\ + d(u, k)) . (A-15) 

This leads to the expression (j20*|l for d(cu, k) and the expression (|2^jl for the 
matrix A. If we equate the imaginary parts, then we discover that 

a(u,k) = e-™{I v - loo), (A-16) 

where 



Neither I n nor are well-defined, but their difference is, and provided one 
expresses a(cu, k) as in (|2*T|). the singularities of these integrals are avoided. 



i rn u)t-2ktan(t/2) 

dt , (A-17) 



dt . (A-18) 



2tt J sin(t/2) 
1 f°° cos(ut - 2k tanh(t/2)) 



2tt J sinh(t/2) 
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